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IX  AMTRACT  (MsjamvmJOOwOftM) 

We  present  in  this  paper  a  study  of  the  azimuthal  anisotropy  of  the  motion  held  observed  during  a  six-day 
campaign  in  March  1986  using  the  MU  radar  in  Shigaraki,  Japan  The  radial  wind  velocity  was  observed  at 
20°  zenith  angle,  at  every  30°  of  azimuth  during  four  days,  and  at  every  43°  dunnz  two  days  A  jet  stream  was 
present  during  the  entire  six  days.  The  average  radial  velocity  variance  from  IU  4  to  19  2  km  was  calculated 
every  four  minutes  and  then  averaged  over  20  minutes  or  one  hour 

The  average  variance  was  found  to  be  a  strong  function  of  both  azimuth  and  time  The  azimuthal  variations 
|  were  analyzed  in  terms  of  the  mean  and  the  first  and  second  harmonics  The  mean  is  proportional  to  the  kinetic 
energy  per  unit  mass  of  the  radial  wind  fluctuations,  and  the  first  harmonic  is  proportional  to  the  vertical  flux 
of  horizontal  momentum  ner  unit  mass.  The  strong  azimuthal  variation  was  usually  dominated  by  the  second 
harmonic;  i.e.,  with  two  peaks,  but  was  occasionally  dominated  by  the  first  harmonic,  with  one  peak  The  phase 
of  the  first  harmonic  was  usually  westward,  but  the  phase  of  the  second  harmonic  was  quite  variable. 

It  was  shown  by  a  development  of  gravity  wave  theory  that  all  of  the  observed  azimuthal  variations  could 
probably  be  caused  by  a  gravity  wave  field  whose  parameters  vary  with  time 
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ABSTRACT 

We  present  in  this  paper  a  study  of  the  azimuthal  anisotropy  of  the  motion  field  observed  during  a  six-day 
campaign  in  March  1986  using  the  Mil  radar  in  Shigaraki.  Japan.  The  radial  wind  velocity  was  observed  at 
20°  zenith  angle,  at  every  20°  of  azimuth  during  four  days,  and  at  every  45°  during  two  days.  A  jet  stream  w  as 
present  during  the  entire  six  days.  The  average  radial  velocity  variance  from  10.4  to  19.2  km  was  calculated 
every  four  minutes  and  then  averaged  over  20  minutes  or  one  hour. 

I  he  average  variance  was  found  to  be  a  strong  function  of  both  azimuth  and  time.  The  a/imuthal  variations 
w  ere  analyzed  in  terms  of  the  mean  and  the  first  and  second  harmonics.  The  mean  is  proportional  to  the  kinetic 
energy  per  unit  mass  of  the  radial  w  ind  fluctuations,  and  the  first  harmonic  is  proportional  to  the  vertical  iluv 
of  horizontal  momentum  per  unit  mass.  The  strong  a/imuthal  variation  was  usually  dominated  by  the  second 
harmonic:  i.e.,  w  ith  two  peaks,  but  was  occasionally  dominated  by  the  first  harmonic,  with  one  peak.  The  phase 
of  the  first  harmonic  was  usually  westward,  but  the  phase  of  the  second  harmonic  was  quite  variable. 

It  was  show  n  by  a  development  of  gravity  wave  theory  that  all  of  the  observed  azimuthal  variations  could 
probably  be  caused  by  a  gravity  wave  field  whose  parameters  vary  with  time. 


1.  Introduction 

Atmospheric  parameters  are  observed  to  fluctuate 
on  all  scales.  When  these  fluctuations  are  approxi¬ 
mately  sinusoidal,  they  can  be  interpreted  as  being  due 
to  waves,  which  in  the  mesoscale  are  almost  certainly 
internal  gravity  waves.  However,  the  fluctuations  are 
sometimes  not  sinusoidal,  so  that  their  cause  is  not 
immediately  evident.  Dewan  (1979)  and  VanZandt 
( 1982)  suggested  that  the  mesoscale  fluctuations  are 
due  to  a  random  superposition  of  gravity  waves,  in 
analogy  with  the  internal  gravity  wave  model  success¬ 
fully  applied  to  very  similar  fluctuations  in  the  ocean 
(Garrett  and  Munk  1972,  1975 ).  Gage  ( 1979),  on  the 
other  hand,  suggested  that  the  mesoscale  fluctuations 
were  due  to  the  energy  spectrum  of  two-dimensional 
turbulence  (a  limiting  case  of  stratified  turbulence). 

The  dynamics  of  these  processes  and  their  dynamical 
effects  on  the  atmosphere  are  quite  different.  In  the 
gravity  wave  field  the  energy  is  thought  to  cascade  to¬ 
ward  small  scales,  where  the  waves  break  down,  with 
at  least  part  of  their  energy  dissipating  into  inertial- 
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range  turbulence.  On  the  other  hand,  in  stratified  tur¬ 
bulence  the  energy  cascades  from  smaller  toward  larger 
scales.  Other  important  differences  are  that  gravity 
waves  propagate  away  from  the  region  of  their  gener¬ 
ation,  so  that  they  may  have  significant  vertical  fluxes 
of  energy  and  momentum,  which,  in  regions  of  wave 
dissipation,  can  cause  large  mean-flow  accelerations 
and  turbulence  intensities  (Vincent  and  Reid  1983; 
Fritts  and  Vincent  1987;  Reid  and  Vincent  1987). 
These  effects  are  now  understood  to  be  important  for 
global  atmospheric  dynamics  and  structure  ( Houghton 
1978;  Lindzen  1981;  Fritts  1984;  Palmer  et  al.  1986). 

It  is  thus  quite  important  to  describe  the  mesoscale 
fluctuations,  to  understand  their  dominant  dynamical 
processes,  and  to  assess  their  effects  on  the  atmosphere. 
This  paper  is  the  first  of  two  papers  in  which  we  ex¬ 
amine  some  of  the  characteristics  of  the  mesoscale  mo¬ 
tion  field  in  the  troposphere  and  lower  stratosphere, 
utilizing  the  unique  capabilities  of  the  MU  radar  in 
Shigaraki,  Japan.  The  MU  radar  can  sound  in  up  to 
16  different  directions  essentially  simultaneously, 
which  permits  the  conduct  of  experiments  that  would 
be  impossible  using  any  other  radar.  We  have  taken 
advantage  of  this  capability  to  measure  the  radial  ve¬ 
locity  as  a  function  of  azimuth  during  a  six-day  cam¬ 
paign  in  March  1986.  These  data  are  used  in  this  paper 
to  study  the  azimuthal  anisotropy  of  the  variance  of 
radial  vclcHty.  Altln  ?’’.muthal  anisr,,*r>rtv  of 
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the  mesoscale  fluctuations  has  been  studied  in  the  me¬ 
sosphere  (Meek  et  al.  1985;  Ebel  et  al.  1987;  Vincent 
and  Fritts  1987;  Manson  and  Meek  1988),  this  is  ap¬ 
parently  the  first  study  in  the  lower  atmosphere.  The 
companion  paper  by  Fritts  et  al.  ( 1989)  studies  the 
vertical  flux  of  horizontal  momentum  and  other  prop¬ 
erties  of  the  motion  field. 

In  section  2  the  experimental  configuration  and  basic 
data  are  described  in  considerable  detail;  only  a  brief 
description  is  given  in  the  second  paper.  The  back¬ 
ground  wind  and  temperature  structure  during  our  ex¬ 
periment  are  presented  in  section  3,  and  the  variations 
of  the  variance  of  the  radial  velocity  versus  azimuth 
and  time  are  described  in  section  4.  In  section  5  a  model 
is  developed  for  the  azimuthal  anisotropy  of  the  vari¬ 
ance  due  to  a  gravity  wave  field,  and  it  is  shown  that 
probably  all  of  the  observed  variations  could  be  ex¬ 
plained  by  a  time-dependent  gravity  wave  field.  Con¬ 
clusions  are  presented  in  section  6. 

2.  Experimental  configuration  and  velocity  data 

The  experiment  was  conducted  using  the  MU  radar, 
located  at  34.85°N,  136. 10°E.  370  m  MSL.  in  Shiga 
prefecture,  Japan.  The  relief  in  this  region  is  strong, 
with  peaks  up  to  800  m  within  10  km  and  up  to  1200 
m  within  30  km.  The  MU  radar  is  a  monostatic  pulsed 
Doppler  radar,  operating  at  a  frequency  of  46.5  MHz 
(wavelength  6.45  m).  The  antenna  is  circular  with  a 
diameter  of  103  m.  resulting  in  a  two-way,  half-power 
beamwidth  of  2. 5°.  Further  details  about  the  MU  radar 
can  be  found  in  Fukao  et  al.  ( 1985). 

The  present  experiment  was  conducted  from  1 122 
JST  (Japan  Standard  Time;  =UTC  +  9  h)  8  March 
1986.  until  1549  JST  14  March  1986  with  only  minor 
interruptions,  for  a  total  of  6  d.  4.5  h.  During  the  first 
four  days,  12  beams  were  directed  at  a  zenith  angle  X 
=  20°  at  azimuths  0  in  steps  of  30°,  and  four  beams 
were  at  X  =  10°  at  0  in  steps  of  90°.  During  the  last 
two  days.  8  beams  were  at  X  =  20°  at  H  in  steps  of  45°. 
and  one  beam  was  directed  at  the  zenith.  However, 
the  studies  in  the  present  three  papers  used  only  the 
beams  at  X  ^  20°.  whose  configuration  is  shown  in 
Fig.  I. 

The  pulse  length  was  300  m  and  the  range  gates 
were  centered  at  ranges  from  7.5  km  to  20.4  km  in 
steps  of  300  m;  at  X  =  20°  the  range  gates  were  centered 
at  heights  above  ground  from  7.05  to  19. 1 7  km  and  at 
altitudes  above  sea  level  from  7.42  to  19.54  km  in  steps 
of  282  m.  The  average  time  resolution  of  the  data  used 
in  our  analyses  was  237  s  (4  min),  with  2256  records. 

3.  Basic  data  and  meteorological  parameters 

Daily  profiles  of  the  zonal  and  meridional  compo¬ 
nents  of  the  L>ackjm>u,'.u  ii  jrizoinai  wind  were  deter¬ 
mined  from  the  iadar  radial  wind  measurements  as 
follows.  The  zonal  component  at  each  altitude  and  time 
was  determined  by  first  averaging  the  radial  wind  mea- 
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Flo.  I.  Configuration  of  the  radar  beams  used  in  this  study.  All 
beams  used  were  al  a  zenith  angle  V  -  20".  The  beams  al  positions 
indieated  by  open  and  closed  circles  were  used  during  the  first  lour 
days  and  the  last  two  days  of  the  experiment,  respectively 

sured  in  the  90°  and  270°  azimuths,  averaging  the  re¬ 
sult  over  each  successive  24-hour  period  centered  at 
2322  JST.  smoothing  the  resulting  profile  with  a  low- 
pass  filter  with  a  cutoff  at  3  km  in  order  to  avoid  con¬ 
tamination  by  waves  with  short  vertical  wavelengths 
and  low  frequencies,  and  finally  multiplying  by  1/ 
sin20°  =  2.92  to  obtain  the  horizontal  component. 
The  meridional  "find  was  calculated  in  the  same  way. 
The  resulting  zonal  and  meridional  wind  profiles,  pos¬ 
itive  toward  the  east  and  north,  are  shown  in  the  left 
and  right  panels  of  Fig.  2a.  with  the  zero  w  ind  for  each 
day  displaced  to  the  right  by  14.6  m  s  '. 

A  zonal  jet  stream  was  present  during  the  entire  six- 
day  period,  with  a  peak  zonal  wind  of  35  ms  1  at 
about  15  km  on  the  first  day,  increasing  to  58  m  s  1 
at  about  1 2.5  km  on  the  third  day.  and  then  decreasing 
to  48  m  s  1  at  10  to  12  km  on  the  sixth  day.  Below 
the  jet  stream  the  changes  were  much  larger  ;  at  10  Km 
the  horizontal  wind  speed  increased  from  about  4  m 
s  1  on  the  first  day.  to  50  ms  1  on  the  third  day.  and 
then  decreased  only  slightly.  As  a  result,  the  jet  stream 
became  much  broader  on  its  lower  side. 

The  background  temperature  profile  in  Fig.  2b  was 
derived  from  loutine  radiosonde  ascents  at  0000  and 
1200  UTC  (0900  and  2100  JST)  from  Shionomisaki. 
100  km  south  of  the  MU  radar.  The  reported  temper¬ 
atures  were  interpolated  to  1000  m  intervals,  averaged 
over  the  entire  six-dav  period,  and  then  smoothed  with 
the  same  low-pass  filter  as  the  w  ind  profiles.  The  back¬ 
ground  profile  of  A2  (where  .V  is  the  buoyancy  fre¬ 
quency).  also  shown  in  Fig.  2b.  was  determined  by 
differencing  the  mean  temperature  profile  over  2000  m. 

The  altitude  of  the  tropopause  determined  from  each 
radiosonde  ascent  increased  from  about  10.5  km  on 
the  first  two  days  to  about  1 2  km  on  the  last  day.  with 
an  average  of  about  1 1  km. 
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Fir; .  2  (a)  Profiles  of  background  zonal  and  meridional  horizontal  wind  averaged  over  24  h 
centered  on  2322  JST  8-12  March  1980.  The  zero  wind  for  each  day  is  displaced  to  the  right  by 
14.6  ms  '.(hi  Profiles  of  temperature  T  and  buoyancy  frequency  squared  N'  averaged  over  the 
entire  six-day  period. 


Figure  3  shows  the  inverse  Richardson  number  Ri  1 
=  S 1  /  N2.  Here  S  is  the  magnitude  of  the  vectoi  slu^.r 
calculated  from  wind  profiles  averaged  over  1  h  and 
differenced  over  300  m  of  range,  so  that  each  box  in 
Fig.  3  is  1  h  long  and  282  m  high,  and  N2  is  the  mean 
over  the  entire  six-day  period  fmm  Fig.  2b.  The  hoxes 
are  coded  as  follows:  black,  Ri-1  >  1;  open,  1  >  Ri-1 
>  >/v.  blank,  '/)  >  Ri  '.  The  zero  time  on  this  and  suc¬ 
ceeding  figures  is  the  beginning  time  of  our  experiment, 
1 122  JST  8  March.  Since  N 2  is  a  mean,  the  variations 
of  Ri  '  with  time  must  be  due  entirely  to  variations 
of  S2.  The  region  of  generally  large  Ri  1  in  the  lower 
stratosphere  between  about  10  and  12  km  was  caused 


by  the  large  shears  on  the  bottom  of  the  jet  stream, 
which  diminished  during  the  six-day  period  (see  Fig. 
2a).  On  the  lop  of  the  j^t  sucam  the  shears  were  much 
smaller  and  N2  was  larger  so  that  Ri  '  was  considerably 
smaller. 

The  richness  and  .  jrir.Oiiity  of  the  i.tdi«ti  veioCuv 
data  measured  by  the  MU  radar  are  illustrated  by  Figs. 
4  and  5,  which  show  time  series  of  low-pass  and  high- 
pass  filtered  data,  respectively.  The  3  db  cutoff  was  set 
at  a  period  of  33  h  to  ensure  that  gravity  wave  motions 
near  the  inertial  period  (20.95  h  at  the  34.85°  latitude 
of  the  MU  radar)  were  not  attenuated  in  the  high-pass 
time  series  and  did  not  contaminate  the  low-pass  series. 
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Fig.  3.  The  inverse  Richardson  number  Ri"1  =  A l2/S2  versus  time  and  altitude.  Here  ,V:  is 
from  Fig.  2b  and  S',  the  vertical  shear  squared,  is  derived  from  the  wind  profiles  measured  by 
the  radar.  Each  box  is  I  h  wide  and  282  m  high.  The  boxes  are  coded  as  follows:  black.  Ri " 1  >  1 : 
open,  1  >  Ri-1  >  blank.  '/3  >  Ri"1. 


In  each  figure,  panels  (a)  and  (b)  present  the  radial 
velocity  measured  by  the  90°  (eastward)  and  0° 
( northward )  beams,  respectively.  The  velocity  scale  on 
the  ordinate  is  for  the  lowest  range  gate  at  an  altitude 
of  7.42  km  and  each  successive  curve  is  for  282  m 
higher  and  is  shifted  upward  by  2  m  s  '  in  Fig.  4  and 
by  3  m  s'1  in  Fig.  5. 

In  the  low-pass  time  series  at  and  below  the  jet  core 
the  zonal  flow  increased  rapidly  between  about  20  and 
60  h  while  the  meridional  flow  underwent  a  large  os¬ 
cillation  from  about  60  to  150  h.  Above  the  jet  core 
both  components  oscillated  with  periods  of  about  two 
days.  The  high-pass  data  below  the  jet  core  also  exhib¬ 
ited  considerable  variability.  Above  the  jet  core  the 
amplitudes  decreased  and  became  less  organized  with 
increasing  height. 

4.  Variance  as  a  function  of  azimuth  and  time 

The  present  study  of  the  anisotropy  of  the  fluctuation 
field  is  based  entirely  on  the  variance  of  the  radial  ve¬ 
locity  as  a  function  of  azimuth  and  time,  which  was 
calculated  as  follows.  At  each  observation  time  and 
altitude  the  radial  wind  profiles  in  the  zonal  and  me¬ 
ridional  planes  were  averaged  over  320  records  =  2 1 .07 
h  ( about  equal  to  the  inertial  period )  centered  on  the 
observation  ume,  and  the  result  was  smoothed  with  a 
low-pass  filter  with  a  cutoff  at  3  km.  Then  the  back- 
giound  prolife  ai  vuvi.  azimuth  was  calculated  trom 
the  zonal  and  meridional  components.  This  profile  was 
then  subtracted  from  the  observed  profile  for  that  azi¬ 
muth  and  time,  and  the  profile  of  variance  was  cal¬ 
culated. 

The  fluctuations  of  the  variance  can  be  reduced  by 
averaging  over  seveial  altitudes  and  times.  In  order  to 
treat  a  reasonably  homogeneous  region,  we  have  av¬ 
eraged  over  only  the  lower  stratosphere;  i.e.,  over  al¬ 


titudes  from  10.43  to  19.17  km  (ranges  from  11.1  to 
20.4  km ). 

Time  series  of  radial  velocity  variance  at  each  azi¬ 
muth  averaged  over  15  records  =  1  h  are  presented  in 
Fig.  6.  The  variance  varied  strongly  versus  time,  with 
significant  peaks  at  about  42.  57.  70.  81.  90.  109.  and 
134  h.  Most  of  the  time  the  azimuthal  variation  was 
dominated  by  the  second  harmonic,  this  was  particu¬ 
larly  evident  during  some  of  the  peaks,  rite  first  har¬ 
monic  was  usuaily  much  smaller  than  the  second,  but 
occasionally  it  was  larger  when  the  second  harmonic 
was  small. 

The  observed  variance  versus  azimuth  in  a  given 
time  interval  can  be  described  succinctly  by  fitting  the 
mean  and  first  two  sinusoidal  harmonics  bv  least 
squares.  This  description  is  also  convenient  for  inter¬ 
pretation.  since  the  terms  can  be  interpreted  physically. 
The  mean  divided  by  sin:20°  is  related  to  the  kinetic 
energy  per  unit  mass  of  the  radial  w  ind  fluctuations. 
According  to  the  analysis  of  Vincent  and  Reid  (  1 983 ). 
the  first  harmonic  is  proportional  to  the  vertical  flux 
of  horizontal  momentum  per  unit  mass,  which  has  an 
important  effect  on  global  atmospheric  dynamics  and 
which  is  treated  in  considerable  detail  in  the  second 
paper  of  this  scries  ( Fritts  et  al.  1990).  All  three  terms 
can  be  interpreted  in  terms  of  gravity  w  ave  theory,  as 
will  be  shown  in  section  5. 

Four  examples  of  the  resulting  variance  versus  azi¬ 
muth  avemved  ev  er  fr.w  records  { ~20  mm  •  U..J  ’*„■ 
harmonic  fits  arc  shown  in  Fig.  7.  Several  features  are 
illustrated.  The  running  mean  azimuthal  curve  cen¬ 
tered  at  record  number  640  was  dominated  by  the  first 
harmonic.  But  at  record  645.  only  20  minutes  later, 
there  was  a  large  second  harmonic  in  addition  to  a 
significant  first  harmonic,  as  ev  idenced  by  the  unequal 
maxima  and  minima.  The  very  rapid  change  from  re¬ 
cords  640  to  645  was  at  the  beginning  of  the  peak  at 
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Flu.  4.  Low-passed  eastward  (a)  and  northward  (b)  time  series  of  radial  velocity  with  the  3  db  cutoff  period  at  33 
h.  The  radial  velocity  scale  is  for  the  lowest  range  gate  at  an  altitude  of  7.42  km  and  each  successive  curve  is  for  282 
m  higher  and  is  shifted  upward  by  2  m  s 


42  h  in  Fig.  6.  Records  660  and  700  were  both  domi¬ 
nated  by  the  second  harmonic,  and  the  first  harmonic 
was  small,  since  the  peaks  were  nearly  equal,  but  the 
ratio  of  the  second  harmonic  to  the  mean  was  much 
larger  in  record  660.  Much  of  this  discussion  will  be 
devoted  to  interpretation  of  the  wide  variety  of  azi¬ 
muthal  behaviors. 

Figure  8  shows  each  successive  fitted  azimuthal  curve 
centered  on  records  440  to  830  (at  times  from  29.0  to 
80.3  h )  with  the  zero  for  each  successive  curve  shifted 
1  (m  s  1 )"  to  the  right.  The  four  examples  in  Fig.  7 
are  indicated  by  arrows  hetween  40  and  50  h  at  the 
margin  of  the  upper  panel.  The  predominance  of  the 
second  harmonic,  particularly  during  the  peaks  at  42, 
57  and  70  h  is  clearly  visible.  Note  also  the  phase 
changes  during  the  42  and  70  h  peaks. 

These  variations  are  presented  quantitatively  in  Fig. 
9.  The  amplitudes  a0  (upper  bold  curve),  a i  (lower 
bold  curve),  and  a2  (thin  curve)  of  the  mean  and  the 


first  and  second  harmonics  are  shown  in  panel  (a). 
The  corresponding  phase  angles  relative  to  north  of  the 
maxima.  (bold)  and  (/>;  (thin),  are  shown  in  panel 
(b).  These  quantities  are  running  averages  over  1  h 
plotted  every  4  min. 

From  12  to  30  h,  varied  with  a  period  of  about  2 
h.  while  <t>  |  changed  by  ±90°  or  1 80°  in  alternate  cycles. 
Similar  periodic  variations  of  a,,  without  significant 
phase  changes,  also  occurred  later,  particularly  from 
80  to  90  h.  After  45  h.  </>,  tended  to  oscillate  about 
270°.  but  there  was  a  tendency  for  «/> ,  to  depart  from 
270°  during  the  maxima  of  to.  There  did  not  appear 
to  be  any  overall  correlation  or  relation  between  the 
phases  </>  i  and  (t>2- 

5.  Comparison  with  gravity  wave  models 

It  was  pointed  out  in  the  Introduction  that  two  dif¬ 
ferent  dynamical  processes  have  been  proposed  to  ac¬ 
count  for  mesoscale  velocity  fluctuations:  gravity  waves 
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Fig.  5.  High-passed  eastward  (a)  and  northward  (b)  time  series  of  radial  velocity  with  the 
3  db  cutoff  period  at  33  h.  The  radial  velocity  scale  is  for  the  lowest  curve  and  successive 
curves  are  displaced  upward  by  3  m  s'. 
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Fig.  6.  Variance  at  each  azimuth  vs  time  averaged  from  10.4  to  19.2  km  and  over  I  h.  The 
lines  of  data  correspond  to  the  beam  positions  shown  in  Fig.  1 . 


and  stratified  turbulence.  Because  their  effects  on  the 
atmosphere  are  quite  different,  it  is  important  to  eval¬ 
uate  the  relative  significance  of  these  processes.  In  this 
section  we  consider  to  what  extent  the  observed  fluc¬ 
tuations  can  be  explained  in  terms  of  gravity  waves. 
We  do  this  by  developing  models  for  the  radial  variance 
as  a  function  of  azimuth  due  to  gravity  waves.  We  find 
that  probably  all  of  the  observed  azimuthal  curves  can 
be  explained  by  a  suitable  field  of  gravity  waves.  Un¬ 
fortunately,  the  theory  of  stratified  turbulence  is  not 
nearly  as  well  developed  as  the  theory  of  gravity  waves, 
so  that  it  is  not  possible  to  show  whether  or  not  the 
observed  variations  could  be  due  to  stratified  turbu¬ 
lence. 


°  0  90  180  270  360  0  90  180  270  360 

AZIMUTH  ANGLE  (°) 

Fic,.  7.  Examples  of  the  variance  versus  azimuth.  The  data  points 
are  averaged  over  five  records  =  20  min.  centered  on  the  labeled 
record  number  and  time.  The  curve  is  the  sum  of  the  mean  and  first 
and  second  harmonics  fitted  to  the  data  points  by  least  squares.  Note 
that  the  panels  labeled  640  and  645  are  centered  only  20  min  apart. 


The  gravity  wave  model  is  developed  as  follows. 
First,  we  develop  a  model  for  the  radial  variance  as  a 
function  of  azimuth  due  to  a  gravity  wave  with  a  given 
frequency,  wavenumber,  and  direction  of  propagation 
that  is  sampled  over  all  phases  of  the  wave.  Even  though 
this  model  treats  only  a  single  wave  and  ignores  the 
effects  of  Doppler  shifting  due  to  the  background  wind, 
the  resulting  model  curves  resemble  the  observed 
curves.  The  model  also  provides  insight  into  how  grav¬ 
ity  waves  can  cause  the  observed  variations.  Second, 
the  effect  of  a  broad  frequency  spectrum  is  explored 
by  integrating  over  a  power-law  spectrum.  Third,  a 
more  general  azimuthal  distribution  of  propagation 
vectors  is  considered.  Finaliy.  we  consider  the  v  ariance 
due  to  quasi-inertial  waves  that  are  not  sampled  over 
all  phases.  Although  the  azimuthal  curves  from  these 
models  resemble  the  observed  curves,  we  do  not  at¬ 
tempt  actually  to  fit  the  model  and  observed  azimuthal 
curves  for  reasons  that  will  become  apparent. 

The  azimuthal  variation  of  the  variance  due  to  a 
gravity  wave  field  can  be  modeled  by  a  slight  general¬ 
ization  of  the  model  that  VanZandt  (  1985 ).  follow  ing 
Pinkel  ( 1981 ).  developed  for  the  power  spectrum  of 
radial  velocity  observed  by  a  Doppler  radar.  The  as¬ 
sumptions  are  the  same  as  those  of  VanZardt  (  1985) 
except  as  noted.  In  particular,  we  have  assumed  that 
there  is  no  background  wind.  For  a  radar  pointing  at 
an  azimuth  0R  =  0°  and  zenith  angle  X.  the  observed 
profile  of  radial  velocity  r,  at  a  range  r  is  given  by 

vt  =  Re  {((-),  sinx  +  <I>  1  cosx  ).  f  ( A  .  «.  0) 

X  exp  [/(At  sinX  cos  fl  ±  mr  cos  X  —  co/ )  ]  |  (  1  ) 

where  the  upper  and  lower  signs  arc  taken  for  upward 
and  downward  phase  propagation  (downward  and  up¬ 
ward  group  propagation ) .  Here  0,  =  cost?  +  /( 1  /w)  sinfl 
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and  <b  1  =  (w/A’){[l  -  (l/ce)2]/[l  ^  («/^)2]},/2. 
where  u >  and  A1  are  the  intrinsic  and  buoyancy  frequency 
normalized  by  the  inertial  frequency,  and  A ( k ,  co,  0)  is 
the  amplitude  of  the  wave.  Variables  k  and  m  are  the 
horizontal  and  vertical  wavenumbers,  and  the  sinx  and 
eosX  terms  are  the  contributions  to  the  radial  velocity 


from  the  horizontal  and  vertical  velocity  components, 
respectively.  Of  course,  ( ! )  is  valid  only  for  1  +  ^  <  V. 

If  we  assume  that  the  gravity  wave  energy  spectrum 
is  variables-separable.  following  Garrett  and  Munk 
( 1975 )  and  Vanz.andt  ( 1982 ).  then  the  spectrum  as  a 
function  of  azimuth  0  is  given  by 


Ex.(m.  cc.  fl)  ~  (A±/2w)EA(m/mt)B( 


co)  J 


i(cos 0  +  /co  1  sinfJ )  sinx  +  <l>  1  cosX 


(2) 


where  A+  and  A  ( =  1  -  A,)  are  the  fractions  of  energy 
associated  with  upward  and  downward  phase  propa¬ 
gation.  respectively;  E  is  the  total  energy  per  unit  mass; 
A(ni/mt)  is  the  wavenumber  spectrum,  where  mn  is 
the  characteristic  vertical  wavenumber:  Z?(w)  is  the 
spectrum  as  a  function  of  normalized  intrinsic  fre¬ 
quency  co:  and  11(H)  is  the  azimuthal  distribution  of 
gravity  wave  propagation  vectors. 

For  the  purpose  of  discussion  we  assume  that  all  of 
the  waves  are  propagating  in  the  azimuth  f)tt  so  that 
//( 0)/2ir  =  5( 0  -  where  b  is  the  Dirac  delta  func¬ 
tion.  Then  the  integration  over  fi  replaces  0  by  f)u-.  But 
since  the  radar  scans  in  azimuth,  it  is  more  convenient 
to  assume  that  the  azimuth  of  the  radar  0H  varies  and 
flu  is  fixed  at  0°.  which  amounts  to  the  transformation 
fl  =  f)lt  -  0H  with  On  =  0°.  Then  with  some  rearrange¬ 
ment 

/•-,.( co.  (IK)  =  EB(ai)  sin:X[5*/2  +  C* con/ 

+  2(1  -2.1  )Cco6.  l/:  cos(lK  +  (S.  /2)  cos20r],  (3) 
where  r>,  =  1  ±  (  l/co):  and  t  =  cotX/A'. 

All  of  the  results  presented  in  this  paper  were  cal¬ 
culated  for  N  =  200.  appropriate  to  the  lower  strato¬ 
sphere  at  the  latitude  of  the  MU  radar.  Then  with  X 


=  20°,  C  =  0.0137.  Equation  (3),  normalized  by  the 
mean,  is  plotted  versus  H  in  Figs.  lOa-c  for  1  =  1. 

3/4,  and  respectively,  and  for  w  =  1.  10 1  :.  10. 
103/2.  I02,  and  200.  (Note  that  in  Fig.  10c  the  curves 
forco=  10l/2and  101/2happer  to  be  almost  identical.) 

The  changes  in  the  relative  magnitudes  of  the  peaks 
can  be  understood  as  follows.  Further  insight  can  be 
found  in.  e.g..  Hines  ( I960.  1974).  Gossard  and  Hooke 
(1975).  and  Pinkel  (  1 98 1  ) .  For  downward  phase 
propagation,  the  gravity  wave  motions  lie  in  a  plane 
of  constant  phase  whose  direction  normal  is  tilted  to¬ 
ward  180°  bv  an  angle  ii  =  tan  '(<l>  ' )  (  =  oc  /  .Y  for  1 
<?  CO  «  A’).  When  U>  =  I.  10 1/2.  10.  !()’  :.  !0;.  or  200. 
fi  =  0°,  0.86°.  2.85°.  9.09°.  30.0°.  or  90°.  respectively. 
Thus  the  phase  plane  remains  nearly  horizontal  until 
co  >  1 .  When  co  =  1  the  motion  is  circular,  so  that  the 
variance  is  the  same  in  all  azimuths.  But  as  co  increases 
the  motion  becomes  nearly  linear  along  the  intersection 
of  the  tilted  phase  plane  and  the  vertical  plane  (ln 
=  0°.  Then  at  co  =  200,  where  (i  =  90°.  the  variance 
is  again  the  same  in  all  azimuths. 

When  the  motion  is  nearly  linear,  the  variance  versus 
t)  is  nearly  proportional  to  cos2i/.  where  4  is  the  dihedral 
angle  between  the  radar  beam  direction  and  the  line 
of  motion  (see  Fig.  11).  At  (I  =  0°  and  180°.  cosy! 
=  cos(7t/2  -  X  ±  fi)  =  sin(x  +  d).  and  at  H  =  90°  or 
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l  ie,  X.  .Successive  titled  curves  centered  20  min  apart  from  2+0  to  80.3  h.  with  the  zero  for  each  successive  curve- 
shifted  to  the  right  by  1  ( m  s  ' )-  The  four  examples  shown  in  Fig.  S  arc  indicated  by  arrows  m  the  margin  of  the 
upper  panel. 
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Tic..  4.  ( a t  1  he  amplitudes  «„  ( upper  hold  curse  ).  ti,  (lower  hold  curve ).  and  n:  ( thin  curve ) 
of  ihe  mean  and  the  first  and  second  harmonies.  ( h)  The  corresponding  phase  angles  e> ,•  ( hold  ) 
and  <.y.  i  thin  )  of  (lie  maxima  relative  to  north.  ./>,  is  not  plotted  when  n,  0.2.  i  hese  quantities 
are  running  avenges  over  I  h  plotted  every  1  mm. 


27()c.  cos i  cos X  sin /i  When  ve  is  small  and  the  mo-  cosij  0°:  i.c..  the  motions  arc  normal  to  the  1  NO ' 
lions  arc  nearly  horizontal.  then  cos\J  -  sinX  at  ft  -  0°  beam.  Thus,  at  larger  w  there  is  only  one  maximum 
and  180°.  which  is  much  larger  than  cosy!  at  H  90°  and  one  minimum,  as  shown  in  T  ig.  10a  for  100. 
and  270".  so  that  the  variance  has  two  maxima  and  I  he  westward  maximum  is  consistent  with  phase  sur- 
t  vvo  minima,  as  shown  in  Fig.  10a  for  ue  -  l()l/:,  10.  faces  whose  normals  are  lilted  toward  the  east  due  either 
anti  1 1)32.  As  w  increases  and  the  line  of  motion  steep-  to  generation  by  shear  instability  below  the  jet  stream, 
cns.  cost/  at  II  -  90".  1 80".  and  270°  decreases  relative  or  to  topographic  forcing  or  convection  at  lower  lex  els. 
to  a  0°.  Indeed,  when  ti  -  X  -  20°  (at  ce  -  68.4),  Of  course,  upward  phase  propagation.  1  0.1.  1. 


■X..  . -  .  L  . . ■-  . 

(a)  A-  =1  (b!  A- =  3/4  (c)  A-  =  i/2 
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Flo  It).  Radial  velocity  variance  versus  radar  beam  a/itnuth  t>  for  a  radar  beam  zenith 
•  ingle  X  20°.  as  a  function  of  gravity  wave  frequency  w,  normalized  by  the  inertial  frequency . 
Ihe  frequency  code  is  m  I.  solid;  l()l,:.  dotted.  10.  dashed;  IO'/!.  dot-dash;  100  and  200. 
‘olid.  T  he  gravity  waves  are  all  assumed  to  be  propagating  in  the  a/imuth  0°.  and  .1  is  the 
traction  of  the  saves  with  downward  phase  propagation  (  upward  group  propagation).  1  Note 
that  in  panel  ( e )  the  curves  for  ie  I ()' and  I01  ’  arc  almost  identical.] 
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I  k..  II  Sketch  illustrating  the  dihedral  angle  between  radar  beams 
and  grants  wase  perturbation  velocities  in  the  vertical  plane  oT II „ 
0  I'm  downward  and  rightward  phase  propagation.  The  double- 
headed  arrows  indicate  the  perturbation  motions  in  three  sur  aces  of 
constant  phase.  The  sketch  is  Tor  radar  beam  zenith  angle  v  r  20° 
and  gravits  wave  frequency  w  51).  resulting  in  an  inclination  angle 
■i  14.5  .  and  dihedral  angles  <1  55.5°  and  H4. 5”. 


produces  the  same  pattern  shifted  by  180°.  so  that  in 
F  ig.  I()c  for  1  =  0.5  the  curve  between  90°  and 

270°  repeats  the  curve  from  270°  to  90°. 

The  frequency  spectra  from  24-hour  time  series  of 
the  present  dataset,  analyzed  in  the  following  paper  by 
F-ritts  et  al.  (  1990).  follow  a  power  law  uc  n  with  p 
*  Therefore,  (3)  has  been  integrated  numerically 
overall  w  from  I  to  200  ( i.e..  from  the  inertial  frequency 
to  the  buoyancy  frequency)  with  p  =  1.  %.  V.,-  and  2. 
The  resulting  curves  are  shown  in  F'ig.  1 2. 

The  curves  in  Figs.  10  and  1 2  are  symmetrical  about 
I) n  .  so  that  a  model  with  all  waves  propagating  in  one 
direction  cannot  explain  the  unequal  minima  that  are 
often  observed  in  Figs.  7  and  8.  In  order  to  demonstrate 
that  a  gravity  wave  model  can  explain  such  minima. 


Fig.  13  shows  the  sum  of  the  curve  in  Fig.  12a  for  p 
=  %  plus  the  same  curve  multiplied  by  '/>  and  shifted 
90°  toward  the  left  (counterclockwise). 

Of  course,  in  reality  the  distribution  of  propagation 
vectors  cannot  consist  of  b  functions.  But  the  fact  that 
the  widths  of  the  maxima  and  minima  in  the  mi  del 
curves  are  similar  to  the  observed  widths  suggests  that 
the  peaks  in  the  actual  11(D)  must  indeed  be  narrow. 

The  foregoing  models  assume  that  the  gravity  waves 
in  the  held  are  sampled  uniformly  over  all  phases,  cither 
versus  space  or  time.  In  the  present  experiment,  av¬ 
eraging  over  the  altitude  range  of  8.74  km  samples  over 
the  phase  of  most  waves,  since  the  characteristic  v  ertical 
wavenumber  in  the  lower  stratosphere  is  only  2-3  km 
( Fritts  et  al.  1988 ).  On  the  other  hand,  averaging  over 
20  min  in  Figs.  7-9  samples  over  the  phase  of  only 
waves  with  short  periods.  For  waves  with  u>  ;•  3,  where 
the  motion  is  approximately  linear,  the  partial  sampling 
affects  the  amplitude,  but  not  the  shape  of  the  curves 
in  Figs.  10,  12.  and  13.  But  for  waves  with  w  *•.  3,  i.e.. 
near  the  inertial  frequency,  according  to  hq.  (  1  )  the 
perturbation  velocity  rotates  on  an  ellipse,  so  that  the 
curves  in  Fig.  10  also  rotate. 

The  rotations  occasionally  seen  in  Fig.  8  may  be  due 
to  such  quasi-inertial  waves  with  long  vertical  wave¬ 
lengths  In  particular,  the  counterclockwise  rotation 
just  after  70  h  could  be  due  to  a  quasi-inertial  wave 
with  a  horizontal  wavelength  of  about  1000  km  that 
has  been  Doppler-shifted  to  a  period  of  about  six  hours 
by  the  background  wind. 

It  is  clear  that  probably  all  of  the  observed  curves  in 
Figs.  7  and  8  and  the  relative  amplitudes  ard  phases 
in  Fig.  9a  and  9b  could  be  explained  by  a  grav  ity  wave 
field  with  suitably  chosen  time-dependent  B(  u. ).  11(D). 
and  .1  .  Indeed,  it  appears  that  most  of  the  observed 
curves  could  be  explained  by  model  curves  w  ith  a  fre¬ 
quency  spectrum  B(  w )  /  a-  '  1  with  p  =  '/i  and  ap¬ 
propriate  //(D)  and  .  I  .  There  are  tw  o  kinds  of  excep¬ 
tions.  F'irst.  when  a,  >  u:.  i.e.,  when  the  observed  vari¬ 
ation  is  roughly  unimodal.  B(u’)  must  be  dominated 
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Flo.  12.  The  curves  in  Fig  10  integrated  over  w  with  frequency  spectra  H( w  1  z  u.-  '  with 
p  1.  solid;  %,  dotted;  '/,  dashed;  and  2,  dot-dash  Other  parameters  arc  as  for  Fig  4 
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Ik,  I  a.  l'hc  sum  of  the  urve  in  I  ig.  12a  Ur/i  %  plus  the  same 
curve  multiplied  h>  ai'vi  shifted  90°  towau  'he  left.  This  process 
i’cnci.  'es  unequal  minima. 

by  relatively  high  frequencies  and  with  .1  =£  '/;•  Sec¬ 
ond.  during  the  periods  of  rotation  /T  ( <  -  I  must  he 
dominated  by  a  quasi-inertial  w-ave. 

We  have  not  attempted  actually  to  In  model  curves 
to  the  observ  ed  curves,  because  the  number  of  param¬ 
eters  is  so  large  that  the  fittings  could  not  be  unique. 
Moreover,  the  present  model  ignores  the  effects  of 
Doppler  shifting  due  to  the  background  wind,  which 
was  quite  large  during  this  experiment.  Doppler  shifting 
has  no  effect  on  the  component  of  motion  at  right  an¬ 
gles  to  the  wind.  Rut  in  the  plane  of  the  wind  it  changes 
the  shape  of  the  frequency  spectrum  ( Sehetiler  and  Liu 
1986;  Fritts  and  VanZandt  1987).  For  radial  spectra 
it  tends  to  reduce  the  effective  value  of/).  Thus,  Doppler 
shifting  would  tend  to  move  the  model  results  toward 
the  curves  with  smaller  />  in  Fig.  12. 

6.  C  onclusions 

We  find  that  the  radial  velocity  variance  averaged 
from  10.4  to  19.2  km  is  a  strong  function  of  azimuth 
Wand  time.  The  azimuthal  variations  were  analyzed  in 
terms  of  the  mean  and  first  and  second  harmonics. 
I  he  amplitude  a  0  of  the  mean  divided  by  sirrx  is  re¬ 
lated  to  the  kinetic  energy  per  unit  mass.  The  first  har¬ 
monic  is  proportional  to  the  vertical  llux  of  horizontal 
momentum,  which  is  discussed  in  the  second  paper 
(  Fritts  et  al.  1 990 )  of  this  series. 

Usually,  the  azimuthal  variation  was  dominated  by 
the  second  harmonic,  but  the  ratio  of  the  amplitude 
of  the  second  harmonic  to  the  mean  and  first  harmonic- 
terms  was  quite  variable.  When  the  first  harmonic  was 
large  enough,  it  often  made  the  two  peaks  and  the  two 
valleys  unequal.  Indeed,  the  first  harmonic  was  occa¬ 
sionally  so  large  that  there  was  only  one  p-’ak.  The 
phase  of  the  maximum  of  the  first  harmonic  was  usu¬ 
ally  westward,  but  the  phase  of  the  second  harmonic 
was  quite  variable. 


It  was  shown  by  a  development  of  gravity  wave  the¬ 
ory  that  probably  all  of  the  observed  azimuthal  vari¬ 
ations  could  be  caused  by  a  held  of  gravity  waves.  I  he 
variations  with  time  must  then  be  due  to  variations  of 
the  parameters  of  the  gravity  wave  held. 

Because  the  theory  of  stratified  turbulence  is  not  suf¬ 
ficiently  well  developed,  it  is  not  possible  to  develop 
model  azimuthal  curves  that  can  be  critically  compared 
with  the  observed  curves.  Thus,  it  cannot  be  deter¬ 
mined  to  what  extent  the  observed  variations  may  be 
due  to  stratihed  turbulence. 

The  present  observations  and  conclusions  strictly 
apply  only  to  the  conditions  under  which  the  obser¬ 
vations  were  made:  the  observations  were  made  in 
March  in  the  middle  latitudes  of  east  Asia,  a  jet  stream 
was  overhead,  and  the  relief  near  the  MU  radar  is 
strong.  In  order  to  understand  how  much  our  results 
depend  on  these  particular  conditions,  similar  exper¬ 
iments  under  other  conditions  are  needed.  The  effect 
of  background  meteorological  conditions  and  season 
can  be  studied  by  the  MU  radar.  The  effect  of  terrain, 
latitude,  etc.,  can  be  explored  by  other  radars  in  dif¬ 
ferent  terrains. 
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